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Abstract 



Cf) ! This is a short review on the thermal, spectral representation in the real-time version 

> 

of the finite temperature quantum field theory. After presenting a clear derivation of 
the spectral representation, we discuss the properties of its spectral function. Two 



applications of this representation are then considered. One is the solution of the Dyson 
equation for the thermal propagator. The other is the formulation of the QCD sum rules 
at finite temperature. 

1 Introduction 

The real- and the imaginary- time versions of the quantum field theory in a medium (i.e., 
at finite temperature and/or density) have somewhat complementary virtues QX|] . The 
real-time version is closer to the conventional (vacuum) field theory, involving no sum 
over frequencies and so requiring no analytic continuation to real energies. However, the 
price to pay for this closeness is the 2x2 matrix structure of all two-point correlation 
functions. Generally speaking, static thermodynamic quantities of a system are usually 
calculated in the imaginary-time version, while the real-time version appears convenient 
for calculating more detailed, particularly time dependent, quantities. 

Here we are concerned with the real-time spectral representation for the thermal, two- 
point correlation functions of local operators. (If this operator is chosen to be just a field 
operator, the representation gives the corresponding thermal propagator.) Although the 
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real-time version is of rather recent development, such thermal spectral representations 
were obtained by Landau @ as early as 1958, after Kallen and Lehmann || derived 
them for the propagators in the vacuum. Earlier Low 0] had derived representations 
for non-vacuum matrix elements in the time component go at fixed space component 
q of the 4-momentum variable g M , conjugate to the coordinate difference of the two 
operators in the matrix element. The Landau representation is of this variety, extending 
such matrix elements to their ensemble average. 

Of course, Landau's derivation of the spectral representation in real-time was techni- 
cally incomplete, as he did not take into account its 2 x 2 matrix structure. The complete 
representation was written by Semenoff and Umezawa []5| in 1983, after Umezawa and 
his collaborators || had established the real-time version. 

In this work we derive in detail the spectral representation for the two-point cor- 
relation functions of local operators, obtaining the symmetry relations satisfied by the 
spectral function. To keep the kinematics simple we consider Lorentz scalar operators. 
As an example we calculate the spectral function to the leading order in scalar field 
theory. We go on to consider two applications of this spectral representation. One is to 
review the well-known reduction of the self-energy matrix to essentially a single function 
by using its factorized structure 0. The other is to formulate the QCD sum rules in 
the real time, finite temperature field theory [Q. 

The derivation of the spectral representation is given in sec. 2, giving an example 
of calculation of the spectral function in sec. 3. In sec. 4 we consider two applications, 
one to Dyson equation and the other to QCD sum rules. Finally in sec. 5 we present a 
summary of the results derived in the paper. 

2 Spectral Representation 

We choose the contour in the complex time plane as originally proposed by Umezawa ||, 
to get a symmetric 2x2 matrix for the free thermal propagator. It consists effectively 
of two segments, one running along the real axis in the positive direction and the other 
parallel to it, but shifted by —i/3/2, in the reverse direction, where (3 is the inverse 
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temperature T. 

Consider the contour-ordered, two-point correlation function of a local operator. 
This operator may, in general, be a composite one, built out of fundamental fields, 
which may be scalar, spinor or vector fields. Thus in QCD the composite operator may 
be any one of the conserved currents, namely the vector current, q{x)^^q{x) or the 
axial vector current, <?( 2)7^75^- <?(£), where q(x) is the doublet field of u and d quarks 
and Tj's are the Pauli matrices. Besides such bosonic operators, one may also have 
fermionic operators like the so-called baryon currents. For the proton this operator is 
e abc [u aT (x)C^^u b (x)}[y tJ ''y 5 d c (x)} D , where a, b, c are the colour indices, C is the charge 
conjugation matrix and D a Dirac index ||. If the operator is not a composite one, 
but just one of the fundamental fields themselves, the correlation function becomes the 
thermal propagator for that field. 

To avoid kinematic complications, we choose in the following a Lorentz scalar opera- 
tor, which may be composite or fundamental and denote it by 0(x). Further we assume 
it to be bosonic. Then the contour-ordered two-point function may be put in the form 
of a 2 x 2 matrix, 

T( r -ri-i( ( T °( x )°(y)) (O(y-0/2)O(x)) \ ( . 

1 V) { {0(x - i(3/2)0(y)) (TO(x - ip/2)0(y - i/3/2)> ) ' 

where T and T denote the usual time and anti-time ordering and (■ ■ ■) denotes the 
ensemble average; thus for the operator O, 

(O) = Tr Oe~ pH /Z, Z = Tr e~ pH , (2.2) 

H being the Hamiltonian of the system and Tr denoting trace over any complete set of 
states. In momentum space, the Fourier transform is denoted by the same symbol, 

T ab (q) = J d 4 ze^ T ab (z), a, b = 1, 2 . (2.3) 

Let us first consider the 11-component of the matrix. We evaluate the trace over 
a complete set of states \m), m = 1,2,.., which are eigenstates of the 4-momentum 
operator P M with eigenvalues (Pm)^. Denoting (p m ) by E m , it becomes 

(TO{x)0{y)) = Z- 1 ^e" /3£; ™(m|TO(a;)0(y)|m) , (2.4) 
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which is a sum over forward amplitudes weighted by the corresponding Boltzmann fac- 
tors. Again inserting the same set of states to extract the co-ordinate dependence of the 
field operators, we get 

(TO(z)O(O)) = Z- 1 e~ pEm (e(z°)e i(pm - pn) - z + e(-z°)e- i(pm ~ Pn > z 



m,n 



x|(m|O(0)|n)| 2 . (2.5) 

It is now simple to work out the Fourier transform ( |2.3| ). The integration over space 
gives rise to 5-functions in 3-momentum, while that over the time variable produces the 
energy denominators. Inserting a 5-function in the energy variables, we may put it in 
the form 

Tn(q) = -L f K (^ML _ ^ML\ , (2 . 6) 

2tt 7-oo \q - q' Q + it q - q' - te J 

where 

M ± (g , q) = Z- 1 e~ pEm (2vr) 4 ,5(g ± p m T p n ) I (m\O(0) \n) | 2 (2.7) 

m,n 

The double sum over states may be converted back to the product of field operators 
giving 

M + (q) = f d 4 ze l(lz (O{x)O{0)) , (2.8) 



and M~(q) having an identical expression with 0(x) and O(0) interchanged. 

In the case of the vacuum expectation value, the two spectral functions, which are 
functions of q 2 only, can be shown to be equal by the use of the causality requirement. 
However, in the present case, where they are functions of go and \q\ (or q 2 and u ■ q in a 
Lorentz covariant framework with w M being the four- velocity of the medium), a similar 
argument does not go through. But we still have two relations connecting M ± , namely 
the Kubo-Martin-Schwinger relation in momentum space, 

M+(g M ) = eP^M-fa) , (2.9) 

and the symmetry relation, 

M + (-q,)=M-(q,) . (2.10) 

These relations are usually obtained from the operator representation (|2.8| ). They may 
also be obtained from the double sum representation ( |2.7| ) : the relation ( 2.10|) is evident, 
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while to get the relation ( |2.9| ) we interchange the dummy indices m, n in any one of 
M ± (q) and use the 5-function to express E m — E n by q . 
Let us now define the spectral function p as 



p(q ,q) = M + (q ,q) - M~(qo, q) = / d 4 ze^([0(x), O(0)]) 



(2.11) 



which, on noting (|2.10|) , is antisymmetric under q^ — > — g M , 



(2.12) 



Also using 



we can express both M in terms of p 



M + {q, 
M-(g p 



- 1 



We now wish to redefine the energy denominators in 
prescription. For this purpose we write it as, 



(2.13) 



with the Feynman ie 



dq'p 
2tt 



— + 7r8(q -q ) 



x {M + (q' ) + M-(q' Q ) - S gn(q' )(M + (q' ) - M~ (<&))} 



(2.14) 



where the q dependence of M ± is suppressed. Here we have written the integrals in 
first as their principal values and then with the indicated ie prescription, the terms with 
^-functions serving to compensate these changes. Folding the range of integration on to 
(0, oo) and using the relations ( |2.12| ) and ( j2.13| ), we get the desired result [], 



T n (q)=if^p(q' ,q) 



+ 2rr5(q' 2 - q 



% 2 - q'o 2 + * e 



l 



(2.15) 



e /%o| _ 1 

The 22-element of the matrix T may be simplified by invoking the translational 



mvanance, 



T 22 (x - y) = 6(x° - y ){O(y)O(x)) + 6(y» - x u )<0(x)0(y)> 



o „0\ 



(2.16) 



*The change of sign of q while using (2.12) is of no consequence here, since it occurs either as q ■ q 
or q ■ p, where p is a 3- vector to be integrated out over its entire range. 
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Repeating the steps similar to above, we get for T 2 2(q) an expression identical to the 
one for Tu(q), except for complex conjugation of its first term in square bracket, 



^22(9) 



d % 



/2 



■p(?o. o) 



2 /2 

Qo -Qo 



'0 2tt 

The 12- and 21- elements turn out to be identical, 



(2.17) 



T 12 (q) = T 21 (q) 



e /3go/2 
' ePi° - 1 



p(qo,q) 



i / dq' 2 p(q' ,q)5(q' 2 - ql) 



oP\qo\/2 



(2.18) 



ePM - 1 

Recognizing the density distribution function and the free propagator in vacuum with 
Feynman boundary condition, 

1 . i 



n 



A 



2 /2 , • 

qo-Qo +^ 



e /%o| _ 1 ' 

in the above expressions, we collect the results for the components as 



%b(qo, q) 



dq'o 



/o 2vr 

where V® b is the free thermal propagator, 



(2.19) 



U ab 



l + n)A + nA* Jn(l + n)(A + A 



n(l + n)(A + A*) nA + (1 + n) A* 



U(\q \)i( ^V(kol), 



(2.20) 
(2.21) 



with 



U 



VT+n y/n 
Ri a/1 + n 



(2.22) 



The matrix U does not depend on the integration variable in ( |2.19| ) and like the free 
propagator, the correlation function also factorizes, 

T 



%b(qo,q) = U(\q \) 



where 



T(q ,q) 



—T 



* U(\q \) 



(2.23) 



dg'o p(qjq) 

2vr q' 2 -ql-ie 



(2.24) 
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Thus the matrix function is given essentially by a single analytic function T(qo, q). This 
integral representation, as such, may not converge, when it would need subtractions. 
However, this point does not concern us in the following. 

As a simple example, let us calculate T(qo,q) for the case where the operator 0{x) 
is a scalar field operator 4>(x), representing particles of mass m. Then from Eq.(2.11) 
we get immediately the spectral function as 

P(Qo, = 2vr sgn(g )%o -q 2 - m 2 ) 

giving 

T(q) = 2 (2.25) 
q l — m z + te 

which is the free propagator function for the scalar field with Feynman boundary con- 
ditions u 

The important point to notice here is that there is only one spectral function p 
giving all the four components of the correlation function and it may be obtained by 
calculating any one of the components, say Tu. But although p = (M + — M~) is (twice) 
the imaginary part of T, it is not so for Tu : As seen from (|2.6|) , (twice) the imaginary 
part of T u is (M + + M~). The two, however, are related through eq. fl2.13| ), 

p{q ,q) = 2tanh(/3g /2) ImT^q^q) (2.26) 

Their real parts are, however, equal, 

ReT(q ,q) = ReT 11 (q ,q) . (2.27) 

Eqs.(2.23 - 24) constitute the finite temperature extension of the Kallan-Lehmann 
representation in vacuum. The breakdown of Lorentz invariance is explicit as it treats 
q as fixed. (As already mentioned, it can be formally restored, but only by introducing 
the four-velocity of the heat bath.) Note the relations (2.26-27) between (the real and 
imaginary parts of) T and the Tu- component. Such relations also exist among the 

components of the self-energy matrix, to be discussed below. 

^The occurrence of —1 instead of i in the numerator of Eq. (2.25) is due to the presence of i in the 
definition of (2.1) of the correlation function. 
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Before leaving this section, we wish to point out that an equation like (2.7) giving the 
spectral function as a double sum over states is not very convenient for its evaluation. 
For, the ^-function in it seems to indicate that even initial and final states with a large 
number m and n of particles would give a contribution in the low energy region provided 
their momentum difference (p m — p n ) is small. Actually most of these particles are from 
the heat bath and do not participate in the interaction with the operators in the two 
point function. Thus it is much simpler to calculate the spectral function from the 
relevant Feynman diagrams, where the thermal propagator properly takes into account 
the particles of the heat bath in the distribution function present in it. Accordingly, we 
turn to the Feynman diagrams in the next section to calculate the spectral function. 

3 An Example 

Let us evaluate the correlation function perturbatively for the case where 

0(x) = <t>i{x)<t>2{x) , (3.1) 

the fields <f>i(x) and <p 2 (x) being two different scalar fields representing particles of masses 
m\ and 777.2- Then the perturbative evaluation for the 11-component follows from 

T u (q) =iJd A x e iqx (Te l J c ^ )d4y Mx)Mx)MQ)MU) ) , (3-2) 

where J0, int ((f)i, (f> 2 ) is the interaction part of the Lagrangian involving (f>i(x) and (f) 2 (x) 
(Fig.l). Usual Feynman rules apply with the vacuum propagators being replaced by the 
thermal ones, 

/ rf 4 xe^(T0,(a:)0 t (O) ) = - - 2 + 2ix5(k 2 - m 2 ) ni (k), % = 1,2 (3.3) 

where rii(k) is the distribution function of particle type i. To leading order the 11- 
component is given by 

V n (q) = -1 [ -^4{(1 + ni)Ax + mA*}{(l + n 2 )A 2 + n 2 A* 2 } (3.4) 

(271") 
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Figure 1: Two point function for 0(x) = <pi(x)4> 2 (x). The dots indicate higher order 
diagrams, whose nature and contribution depend on C int (x) 



where 



1 



Hi ' 

n 2 - 



1 

e f3w 2 _ I ' 

i 



w 2 



k 2 



2 

mf 



(k — q) +777-2 

i 



- A 2 = - -j ^— (3.5) 

fc^ - mf + ?e {k- q) -171% + le 

Each term in the the integrand of ( |3.4| ) has a product of two propagators. Their 
singularities in k are only due to the poles in these propagators. So the integration over 
ko is performed by closing the contour in the k a — plane |, after which the imaginary part 
may easily be read off, 

d 3 k 1 



ImV 



ii 



7T 



(2tt) 3 AwiW 2 



{(1 + 7li)(l + 77 2 ) + 77 1 77 2 }{5(go - Wl ~ W 2 ) + 5(q + Wi + W 2 )} 

+ {(1 + 77i)77 2 + (1 + n 2 )n 1 }{S(q - w x + w 2 ) + S(q + w x - w 2 )} 



(3.6) 



It will be noted that the factors involving the density distribution functions are not 
appropriate for the interpretation in terms of emission and absorption probabilities of 
the particles JlOfl . The desired interpretation follows, if we use the energy conserving 
5-functions to rewrite it as 



ImVu = coth(/5g /2) I 



(3.7) 



where 



71 



d 3 k 



(2tt) 4wiW 2 



(1 + 77x + n 2 ){S(q -w 1 - w 2 ) -5(q + w 1 + w 2 )} 



+ {n 2 - 77i){5(g -Wi + w 2 ) - 5(q +Wt- w 2 )} 



(3. 



*By power counting the integral in (3.2) is divergent. But the divergence resides only in the real 
part. 



10 



Thus from (2.26) and (3.5) we see that the tanh- and coth- factors cancel out in the 
spectral function giving 

p(qo,q) = 2I(q ,q) , (3.9) 

which agrees with the one obtained from the imaginary-time formulation |10[ . Also 
notice that I(q ,q) is antisymmetric under q — > —q in agreement with Eq.(2.12). 

Although we are discussing in this work only the bosonic propagator and bosonic 
intermediate states in its spectral representation, the cancellation of trigonometric func- 
tions is quite general. As long as the propagator is bosonic, there arises a coth- factor 
in the spectral function, even if the intermediate state is a fermion-antifermion system. 
But if we consider a fermionic propagator, the tanh- and coth- factors interchange in the 
expressions analogous to ( |2.26| ) and ( |3.7| ). Thus the trigonometric factors cancel in the 
spectral representation in all cases. Also we note that if there is a chemical potential \i 
in the Boltzmann factor, the argument, f3qo/2 of the trigonometric factors are replaced 
by P(qo — /-0/2, besides other changes. 



4 Applications 

A. Dyson Equation 

The spectral representation is very general in that it is valid for all q 2 . If, however, 
we are interested only in the immediate neighbourhood of the pole in the propagator, 
the Dyson equation is the appropriate tool to work with (Fig. 2). But even here the 
factorizibility of the full propagator, as established by the spectral representation, may 
be used to show that the self-energy matrix appearing in the Dyson equation has actually 
a simple structure, reducing essentially to a single function. 

At finite temperature the Dyson equation for the propagator matrix in momentum 
space is 



V = V - VqYX> , 



(4.1) 
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Figure 2: The series of diagrams summed up in the Dyson equation for the propagator 
where £ is the self-energy matrix. It has the solution 

V- 1 = V 1 + £ (4.2) 



+ UEU (4.3) 



Using the factorized forms for the propagators, it gives 

(l/D \(-(q 2 -m 2 ) 
^ -l/D* J \ q 2 -m 2 j 

It follows that £ must be of the form 

nq) = U- 1 (qo)^ Q _|* Wfo,) (4.4) 

where £(<?) is the self-energy function. Then we have the solution 

D = ^ (4.5) 

q 2 — m 2 — £ + ie 

All the components of the self-energy matrix can now be expressed through the function 
£. From (4.4) we get 

S n = -£* 2 = Ret + i(l + 2n)Imt 



£ 12 = £ 21 = -2i^n(l + n)im£ (4.6) 

These relations should not be interpreted as constraints imposed on the elements of 
the £-matrix by the solution of the Dyson equation. They are, in fact, automatically 
satisfied by the expressions obtained from perturbation theory. The above method of 
solution is just a convenient way to establish them. 

Let us note here that the factorization ( |4.3| ) of the full matrix propagator is really 
not necessary to solve the Dyson equation. Indeed, writing £ = U^H'U^ 1 and using 
the factorization of the free propagator alone, eqn. ( |4.2| ) gives 
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On obtaining the inverse, one finds that each of the components of the propagator matrix 
is a linear combination of the pole terms, {q 2 — m 2 — E + ) 1 and (q 2 — m 2 — E_) \ where 

S± = - S 22 ) ± i((S' n + E 22 ) 2 - 4S' 12 S 21 ) 1/2 . (4.8) 

Thus although we can solve Dyson's equation, we cannot relate the components of E to 
a single function E, without using factorizibility of the propagator. But on using these 
relations, we see that E± indeed reduces respectively to E and E*. 

The two-particle contribution to the spectral function considered in sec. 3 may also 
be taken as a contribution to the self-energy. Thus the calculation of D u there also 
serves as an example of En. 



B. Thermal QCD Sum Rules 

A vacuum QCD sum rule may be derived for a two-point correlation function, 
(0\TOi(x)O 2 (0)\0) of any two operators 0\(x) and 2 (x), which can be built out of 
quark and gluon fields of the QCD theory. As already mentioned, these operators can 
be mesonic, like the conserved currents of the flavour SU(2) group, or can be fermionic, 



like the baryon currents. The QCD sum rules [JTTJ] are obtained by equating the Fourier 
transforms of the spectral representation to that of the Operator Product Expansion 
(OPE) of the two-point function at a convenient space-like momentum. 

These vacuum sum rules have already been extended to the medium at finite tem- 
perature (with or without chemical potential) 0, 0. It is not our purpose here to 
review this topic in detail. Instead, we merely show how one deals with the 2x2 matrix 
structure of the two-point function in the formulation of these sum rules. 

To bring out the essential features of these thermal QCD sum rules without getting 
involved in kinematic complications, we again consider the two-point correlation function 
(2.1) of a bosonic Lorentz scalar operator 0(x), satisfying the spectral representation 
(2.23). Further, if we take q 2 space-like, q 2 = q\ — q 2 = —Q 2 , Q 2 > 0, where the function 
T(Q 2 ,q 2 ) is real, it reduces to 

T ab {Q 2 ,q 2 )=T{Q 2 ,q 2 ) ( J M (4.9) 
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where 

J-q 2 2,71 q z + Q z 

with p given by (2.26). But as we already pointed out at the end of sec. 2, it is more 
useful to calculate it perturbatively within the context of an appropriate effective field 
theory. 

It is now clear that the 11- (or 22-) component of the matrix (4.9) contains all the 
information in the space-like region. We may then obtain the thermal spectral sum rules 
by considering this component alone. 

The other ingredient needed to write the sum rules is the OPE for the 11-component. 
It is in deriving this expansion that the actual quark and/or gluon structure of the 
operator 0(x) must be spelt out. What the OPE does is to expand the two-point function 
in terms of local operators having finite matrix elements, its singular behaviour at short 
distance being segregated in c number (Wilson) coefficient functions.^ At sufficiently 
high Q 2 , its Fourier transform may be written as, 

(Pi) 

(Q 2 

where c$ are numbers. The local operators Oj are labelled by the index i denoting 
their dimensions. The inverse powers of Q 2 are determined by dimensional counting, 
m = i/2 + 2 — p, where p is the dimension of the operator 0(x). 

The sum rule follows from equating the spectral representation (4.10) to the OPE 
(4.11) in a convenient region of Q 2 . After Borel improvement fil l, it finally becomes 

^/;^-/"V(,fl=E, (m . ( g , ( > Mi) . , (4.i2) 

where M is the Borel mass replacing Q. 

In writing sum rules of practical interest, we have, of course, to deal with more 
complicated kinematics depending on the nature of the operators 0(x) in the two-point 
function. But it is clear that the method of derivation and the structure of these sum 



i J d*xe^(TO(x)O(0)) n = E^tSw , (4.11) 



§ The off-diagonal elements are actually regular as the two arguments of the operator O cannot tend 
to each other for finite /3. 
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rules will be the same as in the above illustrative example. The special features of these 
sum rules over the corresponding vacuum ones and their evaluations are described in 



the literature [12, pL H. 



5 Conclusion 

We present a detailed derivation of the spectral representation and discuss properties 
of the spectral function. We calculate the spectral function, as a proto-type example, 
in a scalar field theory and note its agreement with the expression obtained from the 
imaginary time version of the finite temperature field theory. Relationships similar to 
those among the components of the self-energies are shown to exist also among the 
components of the two-point correlation matrix. 

We then review two important applications of this representation. One is the solu- 
tion of the Dyson equation for the thermal propagator. The factorizability of the full 
propagator (as well as of the free one) leads to a simple solution, relating at the same 
time all the four components of the self-energy matrix to a single function. We also point 
out that the Dyson equation could be solved even without this input of factorizability, 
but then the simple form of the self-energy matrix could not be exposed. The latter 
could, of course, be inferred from its perturbative calculation in each order. 

The other application relates to the formulation of the QCD sum rules in the real 
time thermal field theory. We show how the matrix structure of the two-point function 
simplifies at space-like momenta, where it is needed for writing the sum rules. We 
can then work with its 11- component, obtaining the familiar sum rules written in the 
literature. 
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